Polarized Fermi condensates with unequal masses: Tuning the tricritical point 
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We consider a two-component atomic Fermi gas within a mean-field, single-channel model, where 
both the mass and population of each component are unequal. We show that the tricritical point at 
zero temperature evolves smoothly from the BEC- to BCS-side of the resonance as a function of mass 
ratio r. We find that the interior gap state proposed by Liu and Wilczek is always unstable to phase 
separation, while the breached pair state with one Fermi surface for the excess fermions exhibits 
differences in its DoSs and pair correlation functions depending on which side of the resonance it 
lies. Finally, we show that, when r > 3.95, the finite temperature phase diagram of trapped gases 
at unitarity becomes topologically distinct from the equal mass system. 



Recent advances in the ability to manipulate and con- 
trol ultra-cold atomic vapors have provided a unique ex- 
perimental system in which to explore pairing phenom- 
ena. Following the successful realization of the crossover 
from the BCS to the Bose- Einstein condensed phase [![, 
attention has turned to the consideration of more exotic 
Fermi superfluids. A subject that has attracted partic- 
ular interest is that of Fermi condensates with imbal- 
anced spin populations 0, 0], having relevance to QCD 
and magnetized superconductors [4j. Equally intrigu- 
ing is the case in which both the mass and population 
of each fermionic species in the condensate are unequal. 
Indeed, the realization of Feshbach resonances in Fermi- 
Bose mixtures @ , and the predicted stability of diatomic 
molecules close to resonance over a wide range of mass 
ratios 6], suggest that such superfluid mixtures should 
be experimentally accessible. 

Previous studies 0, @, @, E3] have raised several impor- 
tant issues unique to Fermi condensates with spin and 
mass imbalances. Firstly, it has been proposed that the 
breached pair (BP) state, where the superfluid and excess 
fermionic states phase separate in momentum space, can 
possess excess fermions with two Fermi surfaces (BP-2) 
— the interior gap state Q . In the case of equal masses, 
the BP state can only have one Fermi surface (BP-1). 
However, it is still unclear whether BP-2 can become sta- 
ble for large mass ratios near the Feshbach resonance [ll[ , 
or whether it is always unstable to phase separation in 
real space as in the weak coupling limit [T3| • Secondly, al- 
though there have been studies of the zero temperature 
phase diagram for the homogeneous gas with unequal 
masses so far no tricritical point (such as that dis- 

cussed in the equal mass system F^, 14 1 ) has been identi- 
fied. It is natural to ask whether such a tricritical point is 
generic and, if so, how it evolves with mass ratio. Finally, 
it has been shown that trapped Fermi gases with unequal 
masses can exhibit spatial phase separation at zero tem- 
perature that differs qualitatively from that of the equal 
mass case 10]. Whether and how such features extend 
to finite temperatures remains unanswered. Focussing 
on the homogeneous system, all three issues will be ad- 



dressed in this work. To focus our discussion, we will 
address the phase boundaries between spatially homo- 
geneous phases, leaving the potential for Fulde-Ferrell- 
Larkin-Ovchinnikov phases [l5| to future investigation. 

Referring to wide (viz. entrance-channel dominated) 
Feshbach resonances [3], we restrict attention to a 
single-channel Hamiltonian with contact potential: 

o-=T4 ko- 

+ y X! C k+q/2T C -k+q/2j, C -k'+q/2i C k'+q/2T . (1) 
k,k',q 

Here, eko- = 5^— (ft = 1), g is the interaction strength, V 
is the volume, and we allow the mass m a and chemical 
potential [i a of each spin to be different, with average 
chemical potential fi = (/if + /ijJ/2 and 'Zeeman' field 
h = (/if — Atj.)/2. To obtain the topology of the phase di- 
agram spanning the BCS-BEC limits for unequal masses 
(including the locus of the finite-temperature tricritical 
point at which the transition between the superfluid and 
normal phases switches from first- to second-order), we 
will develop a mean-field analysis of the system analo- 
gous to that presented for equal masses [l3j]. While such 
a treatment is not expected to be quantitatively correct 
for all interaction strengths, it should provide a reliable 
qualitative description even close to unitarity. 

Setting £k+ = t~ — , — = — — I — — , ek- = £k+ 'f , ll , 
and r — the free energy density can be expressed as 
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where minA gives the global minimum with respect to A 
(to be calculated numerically). Here, £k = £k+ — M and, 
defining = y^-j- A 2 , the quasiparticle energies are 
given by E\^ a = qp (h — ek-)- The free energy ([2|) 
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FIG. 1: (Color online) Zero temperature phase diagrams for 
different mass ratios r = my/m^, where a positive polariza- 
tion m corresponds to an excess of particles with spin f. All 
mass ratios exhibit the same three phases: the normal state 
(N), the spatially homogeneous superfluid (SFm), and phase 
separation (PS) between normal and superfluid states in real 
space. The region of SFm expands with increasing r. The 
boundaries enclosing the PS region are all first-order, while 
the tricritical point is represented by a circle. The line defined 
by m/n = corresponds to the usual BCS-BEC crossover. 



differs from that of the equal mass system only through 
the appearance of a momentum-dependent contribution, 
£k-! to the Zeeman term which leads to the symmetry 
(h,r) I— ► (— h, 1/r). We introduce the s-wave scatter- 
ing length a via the prescription = + X)k° 2e^~> 
where the UV cutoff fc can be sent to infinity at the end 
of the calculation. With experiments performed at fixed 
density, fi and h are determined from the total density 
n = n-f + ni = —dOP/d^i and the population imbalance 
m = n-f — n± = —dQ°/dh [22]. Note that the stability cri- 
terion adopted in Refs. [8, 9] corresponds only to finding 
local minima with respect to A in Eq. {2} and thus will 
not correctly locate a first order transition. The presence 
of such a transition implies a region of phase separation 
in real space if n and m are held fixed The two 

phases in question are a normal phase with A = and a 
superfluid phase with A^O, which is less magnetized as 
pairing tends to enforce equality of populations. 

Beginning with zero temperature (Fig.[T|), we find that 
the basic structure of the phase diagram (^r^, ^) for 
unequal masses mirrors that of the equal mass system 
(23| . Phase separation (PS) is found between the normal 
(N) and superfluid states. On the BCS side the super- 
fluid component is always unpolarized, while the BP or 
magnetized superfluid (SFm) state exists for sufficiently 
strong interaction and eventually undergoes a second- 
order transition to the fully-polarized normal state be- 
yond the tricritical point. However, the position of the 
tricritical point is dramatically shifted, with the regions 
of PS and SFm shrinking and expanding, respectively, as 




FIG. 2: (Color online) Evolution of the tricritical point 
(1/kwa, h/ep)tiic as a function of mass ratio r with h > 0. 
The values of r = 1, the equal mass case, and r = 3, where 
the chemical potential fi changes sign, are explicitly marked. 
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FIG. 3: (Color online) Quasiparticle dispersions £kf, Ski 
together with the majority species (|) projection of the 
DoS(oj) = — i Q ! [(?ri(k,ien = u> + irj)] for mass ra- 
tios r = 0.5 (upper panels) {A/\fi\ — 0.25, h/\fi\ = f.06, 
1/kFa = 2.93, m/n = 0.f2) and r = fO.O (lower panels) 
(A/n = 0.85, h/\n\ = f .44, l/k F a = -0.07, m/n = 0.f6). 



r increases (in qualitative agreement with Refs. [8l 15, [icj| ) . 

We can gain further insight into the phase diagram 
by examining the behavior of the tricritical point as a 
function of mass ratio, r. As shown in Fig. [21 the tri- 
critical point, determined as the point at which both the 
quadratic and quartic terms in the Landau expansion of 
the free energy (F2J) vanish, evolves smoothly from the 
BEC to BCS limits, with (h/eF, l/fepcOtric ~~ > i 00 ! 00 ) a $ 
r — > 0, and (h/eF, l/hpa)^ — > (2 2 / 3 ,— oo) as r — * oo. 
Moreover, as in the equal mass system, the tricriti- 
cal point always corresponds to a fully-polarized state, 
m/n = 1 [24|]. It is interesting to note (Fig. ^ that the 
chemical potential, /1, at the tricritical point becomes 
positive for r > 3, hinting at the possibility of a BP-2 
state. However, a mean-field analysis of the phase bound- 
aries in the limit r — » 00 shows that only the BP-1 state 
is ever stable. As in the weak-coupling limit fl2T | , ph ase 
separation is always favored over the BP-2 state [25j |. 
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FIG. 4: (Color online) Averaged momentum distribution for 
each spin state rikf.J. an d correlation function Cfj. for r = 
0.5 (left panels) and r — 10.0 (right panels) and the same 
parameter values used in Fig. [3] 

Even though the BP state only ever possesses one 
Fermi surface, the associated quasiparticle excitation 
spectrum can display minima at non-zero k, a property 
which is also characteristic of the BP-2 phase. As demon- 
strated in Fig. [3l while the spectrum remains gapless for 
both r = 0.5 and r = 10 (i.e., the density of states, 
DoS(cj), is finite at u> = 0), the excitation spectrum for 
r = 10 also exhibits stationary points at non-zero k re- 
sulting in a square-root singularity in the associated DoS. 
Such a singularity may be regarded as 'BCS-like', while 
its absence for r = 0.5 is 'BEC-like' (cf. quasiparticle ex- 
citation spectra at the BCS-BEC crossover in the equal 
mass case |l7j)- More generally, DoS singularities in the 
BP state only occur for r > 1, where the species with 
the smaller mass is in the majority. When /x > 0, the 
singularities exist for all polarizations, but when \x < 0, 
they are restricted to the region around m = where 
A/|/i| > 2y / r/(r — 1) [3]. Such dramatic differences in 
the DoS should be accessible experimentally using optical 
excitations [l~8l |. 

Further signatures of breached pairing are visible in the 
momentum distribution, n k(T , and correlation function, 

C n (ki,k 2 ) = (n klT ra k2i ) - (n klT )(n k2i ) 

= 5 kl ,-k 2 ^[i - e(-£ kT ) - e(-£ ki )] 2 . 

Referring to Fig.|4l breached pairing is characterized by a 
phase separation in momentum space between the excess 
of majority species | and the minority species J, involved 
in the superfluid state. In both cases, the correlation 
function C-j-j, shows a 'hole' for momenta less than the 
Fermi momentum of the majority quasiparticlcs. This 
hole at small momenta is reminiscent of the Pauli block- 
ing observed in the closed-channel molecule of the 40 K 
Feshbach resonance fl9| where there is always an inherent 
particle number asymmetry in the open channel, even in 
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FIG. 5: (Color online) Phase diagrams in the (T/fi,h/^j,)- 
plane for different mass ratios r at unitarity 1/kpa = 0. Mass 
ratio r = 6.7 corresponds to the case of 6 Li- 40 K mixtures. The 
transition between the N and SF phase can be first (red) or 
second order (black solid), separated by the tricritical point. 
The dashed black line describes the second order transition 
in the region where the true transition is first order. Lines of 
constant gradient T/h represent the trajectory traveled from 
the center to the edges of the trap, examples of which are 
shown by the arrows (a-e). The starting points on these lines 
are the values of (T/fi, hj \j) at the trap center. 

the usual BCS-BEC crossover. The correlation function 
C\\ can also directly probe the sign of //: for \i > we 
have a peak beyond the blocked region, as shown in the 
r = 10.0 case, provided A/fi > \h/pb+(l-r)/(l+r)\. In 
principle, such a feature can be detected experimentally 
using noise correlations (20j . 

Turning to finite temperature, the phase diagram for 
unequal masses in the homogenous system is qualitatively 
similar to the equal mass case, but the topology of phases 
within a trap can be different. A study of 6 Li- 40 K mix- 
tures at zero temperature [l(| has already revealed that 
the unequal mass case offers a richer variety of phase- 
separated states. In addition to configurations where 
there is a superfluid core surrounded by the normal state 
(as in the equal mass case), the superfluid region can oc- 
cupy a shell sandwiched between a normal inner core and 
a normal outer shell with opposite polarizations. To in- 
vestigate the effects of temperature T and different mass 
ratios r in the trap geometry, we restrict ourselves to 
the case in which the local density approximation can 
be applied. Thus, assuming each species experiences the 
same trapping potential V(r), the effects of the trap can 
be included in the spatially-varying chemical potential, 
fi — V(r), while the term h remains fixed. 

In Fig. [5l we plot mean-field phase diagrams at uni- 
tarity as a function of T//J, and h/fi. Note that Nozieres- 
Schmitt-Rink fluctuations will not alter these diagrams 
since corrections only enter into the density, n, and the 
polarization, m [2l[. We restrict ourselves to the case 
where /x > because this is enough to completely en- 
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compass the SF-N phase boundaries, even though a fully- 
polarized normal state obviously exists for —\h\ < /i < 0. 
While the equal mass (r = 1) phase diagram is symmet- 
ric around ft = as expected, we see that lowering or 
increasing r translates the superfluid dome to the left or 
right, respectively. For sufficiently small translations, we 
have the ordinary trapped case of a superfluid core with 
a first (arrow a) or second (arrow b) order transition to 
the surrounding normal state. However, once r > 3.95, 
the superfluid region has shifted entirely to the h > 
plane, as shown for r = 6.7 and r = 10 27]. Provided 
h//j,is sufficiently small at the trap center, this naturally 
leads to a superfluid shell structure where the superfluid 
phase is sandwiched between a 'heavy' normal core and 
an outer 'light' normal phase. This structure can either 
have two first-order SF-N phase boundaries (d) , one first 
and one second (c), or two second-order phase bound- 
aries (e), depending on the value of T/h. The case con- 
sidered in Ref. (Toj corresponds to the T/fi = axis in 
the r = 6.7 phase diagram, and is thus an example of 
category d. Note, however, that one tricritical point van- 
ishes when r is further increased (corresponding to the 
point where (l/&Fo)tric changes sign in Fig. [2j, so that 
'd-type' SF shells are then no longer possible. 

In conclusion, we have investigated the mean-field 
phase diagram of a polarized Fermi condensate with un- 
equal masses. We have shown that the zero temperature 
tricritical point smoothly crosses over from 1/kp-a > to 
l/kpa < as r increases, but the interior gap state is 
never stable, even in the limit of infinite r. However, dif- 
ferences in the BP states do show up in the DoS and pair 
correlations. Finally, we show how the phase diagram of 
trapped gases depends on r and T, including how one 
obtains superfluid shells for r > 3.95. 



* Electronic address: mmp24@cam.ac.uk 

* Electronic address: fmm25@cam.ac.uk 

[1] C. A. Regal et al, Phys. Rev. Lett. 92, 040403 (2004); 
M. W. Zwierlein et al, Phys. Rev. Lett. 92, 120403 
(2004); C. Chin et al, Science 305, 1128 (2004); T. Bour- 
del et al, Phys. Rev. Lett. 93, 050401 (2004); J. Kinast 
et al, Phys. Rev. Lett. 92, 150402 (2004); M. Zwierlein 
et al, Nature 435, 1047 (2005). 

See, e.g., P. F. Bedaque et al, Phys. Rev. Lett. 91, 
247002 (2003); J. Carlson and S. Reddy, Phys. Rev. Lett. 
95, 060401 (2005); C.-H. Pao et al, Phys. Rev. B 73, 
132506 (2006); D. T. Son and M. A. Stephanov, Phys. 
Rev. A 74, 013614 (2006); T. Mizushima et al, Phys. 



[2 



Rev. Lett. 94, 060404 (2005); D. E. Sheehy and L. Radz- 
ihovsky, Phys. Rev. Lett. 96, 060401 (2006). 

[3] M. W. Zwierlein et al, Science 311, 492 (2006); G. B. 
Partridge et al, Science 311, 503 (2006); M. W. Zwierlein 
et al, W. Ketterle, Nature 442, 54 (2006; Y. Shin et al, 
Phys. Rev. Lett. 97, 030401 (2006). 

[4] R. Casalbuoni and G. Nardulli, Rev. Mod. Phys. 76, 263 

(2004) . 

[5] C. A. Stan et al, Phys. Rev. Lett. 93, 143001 (2004); 

S. Inouye et al, Phys. Rev. Lett. 93, 183201 (2004). 
[6] D. Petrov et al, J. Phys. B: At. Mol. Opt. Phys. 38, S645 

(2005) . 

[7] W. V. Liu and F. Wilczek, Phys. Rev. Lett. 90, 047002 
(2003). 

[8] M. Iskin and C. A. R. Sa de Melo, Phys. Rev. Lett. 97, 

100404 (2006); cond- mat/0606624. 
[9] S.-T. Wu et al, Phys. Rev. B 74, 224504 (2006). 
[10] G.-D. Lin et al, Phys. Rev. A 74, 031604(R) (2006). 
[11] M. M. Forbes et al, Phys. Rev. Lett. 94, 017001 (2005). 
[12] P. F. Bedaque et al, Phys. Rev. Lett. 91, 247002 (2003). 
[13] M. M. Parish et al, Nature Physics 3, 124 (2007). 
[14] K. B. Gubbels et al, Phys. Rev. Lett. 97, 210402 (2006). 
[15] P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 (1965); 
A. I. Larkin and Yu. N. Ovchinnikov, Sov. Phys. JETP 
20, 762 (1965). 

[16] T. Kohler et al, Rev. Mod. Phys. 78, 1311-1361 (2006). 
[17] M. M. Parish et al, Phys. Rev. B 71, 064513 (2005). 
[18] P. Torma and P. Zoller, Phys. Rev. Lett. 85, 487 (2000). 
[19] M. M. Parish et al, Phys. Rev. Lett. 94, 240402 (2005); 

J. M. Acton et al, Phys. Rev. A 71, 063606 (2005). 
[20] E. Altman et al, Phys. Rev. A 70, 013603 (2004). 
[21] P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 

195 (1985). 

[22] The definition of m also leads to the symmetry (m, r) i— > 
(— m, 1/r), allowing restriction to m > when varying r. 

[23] We define the Fermi momentum with respect to average 
density n/2, so that &f = (37r 2 n) 1//3 and ef = k F /2m r . 

[24] This feature survives Nozieres-Schmitt-Rink fluctuation 
corrections, as in the equal mass case but other 

corrections that include interactions between excess 
fermions and molecular bosons may modify it. 

[25] The existence of a BP-2 phase for r < 1 was discussed 
in Ref. 0], but the possibility of phase separation was 
ignored. 

[26] For zero polarization, while it is true that the ground 
state remains insensitive to mass ratio (as emphasized in 
Ref. [9J), the quasiparticle excitation spectrum and its as- 
sociated DoS are strongly dependent on it. In particular, 
only one branch of the spectrum has a minimum at non- 
zero k when |/x/A| < \r — 1|/2a/t, while both branches 
have minima at non-zero k when |r — l\/2^/r < [i/ A and 
neither branch possesses them when fi/A < —\r—l\/2y/r. 

[27] Equivalently, the whole superfluid region will shift to the 
h< plane if 1/r > 3.95. 



